Modeling spin transport in electrostatically-gated lateral- channel silicon devices: role 

of interfacial spin relaxation 
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Using a two-dimensional finite-differences scheme to model spin transport in silicon devices with 
lateral geometry, we simulate the effects of spin relaxation at interfacial boundaries, i.e. the exposed 
top surface and at an electrostatically-controlled backgate with Si02 dielectric. These gate- voltage- 
dependent simulations are compared to previous experimental results and show that strong spin 
relaxation due to extrinsic effects yield an Si/Si02 interfacial spin lifetime of « Ins, orders of mag- 
nitude lower than lifetimes in the bulk Si, whereas relaxation at the top surface plays no substantial 
role. Hall effect measurements on ballistically injected electrons gated in the transport channel yield 
the carrier mobility directly and suggest that this reduction in spin lifetime is only partially due to 
enhanced interfacial momentum scattering which induces random spin flips as in the Elliott effect. 
Therefore, other extrinsic mechanisms such as those caused by paramagnetic defects should also be 
considered in order to explain the dramatic enhancement in spin relaxation at the gate interface 
over bulk values. 
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In nonmagnetic materials, a nonequilibrium spin po- 
larization will relax away over a characteristic timescale 
called the "spin lifetime" . If the material has weak spin- 
orbit interaction, dilute nuclear spin species, and lattice 
spatial inversion symmetry, the conduction electron spin 
is relatively insulated from other angular momentum de- 
grees of freedom and the relaxation rate is low, resulting 
in long spin lifetime Sihcon (Si) has these prop- 

erties and the predicted long conduction electron spin 
lifetime in the bulk has been confirmed with magnetic 
resonance linewidth measurements Q and in recent trans- 
port experiments probing spin time-of-flight with spin 
precession [5, 6]. 

This nearly spin-conserving property of bulk Si is 
encouraging for the development of applications such 
as schemes to utilize spin to encode and process 
information. However, the spin-enabled devices and 

circuits necessary in any implementation will invariably 
require surfaces and interfaces where the circumstances 
preserving spin in the bulk no longer apply. For instance, 
paramagnetic defects, 12 1 inversion asymmetry, ^is'] and 
fluctuating electric fields, IJ] all of which are absent in 
the ideal bulk, are present at interfaces and can drasti- 
cally increase the spin relaxation rate resulting in a much 
shorter spin lifetime. [iMl 



In prior work, we have experimentally shown a sub- 
stantial spin polarization suppression controlled by the 
attractive voltage bias on an electrostatically gated 
Si/Si02 interface in long-distance lateral- geometry de- 
vices making use of ballistic hot electron transport 
through metallic thin film contacts for spin injection 
and detection. Q Simultaneously, we observed a reduc- 
tion in the mean and standard deviation of the spin 
transport time obtained through coherent spin precession 
measurements. [6] This behavior is in contrast to vertical 
spin transport devices under similar conditions where a 
shorter transit time always results in a larger spin polar- 



ization because of the finite bulk lifetime. Q 

These counterintuitive results were analytically repro- 
duced with a one-dimensional model in Ref. |3 by simu- 
lating the gate voltage dependence of both spin polar- 
ization and spin transport time measurements with a 
phenomenological effective spin lifetime which decreased 
with increasing gate voltage. Although the experimen- 
tal data could be matched successfully with the results 
of this simple scheme, the detailed role of the Si/Si02 
interface (and exposed air/Si top surface) in spin relax- 
ation and its intrinsic spin lifetime distinct from the bulk 
could not be discerned. To accomplish this, it is neces- 
sary to use a two dimensional modelfl9l] since the gate 
is expected to change the transverse proximity of trans- 
ported spin-polarized electrons to the Si/Si02 interface, a 
detail which cannot be captured by any one-dimensional 
model. 

In the present work, we develop a two-dimensional 
finite-differences model to simulate transport in electro- 
statically gated lateral spin transport devices. Hall mea- 
surements on ballistically injected electrons are used to 
provide input on the effects of interfacial momentum 
scattering affecting charge transport. We show that the 
strong spin polarization suppression accompanying a re- 
duced mean transit time and uncertainty observed in 
spin transport experiments is indeed accounted for by 
a Si/Si02 interfacial spin lifetime of ~ Ins, smaller than 
the bulk Si value by at least two orders of magnitude. 0] 
Furthermore, inclusion of spin relaxation at the exposed 
air/Si surface is shown to produce incompatible trends in 
the spin transport characteristics and is therefore inade- 
quate to account for the dominant mechanisms responsi- 
ble for our empirical observations. 

Our model geometry is shown in Fig. 1(a). The lat- 
eral transport channel is lO/xm thick, and the spin in- 
jector and detector contacts are 10/im long. The total 
lateral extent of our computational region is 100 /im. 
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FIG. 1: (a) Finite-difference computational domain showing 
location of spin injector, detector, and gate, (b) Example 
electrostatic potential landscape (solution of Poisson Eqn. JT] 
plus S.leV conduction band offset between Si and Si02),[3| 
through which we model spin-polarized transport according 
to Eqn. [2l Here, Vg=2V and yD=10V relative to the injec- 
tor potential. Spin relaxation is nonzero only at either the 
top or bottom interfaces as explained in the text; insulating 
boundary conditions are used except at the detector where 
spin flows out to consititute a spin current, matching experi- 
mental conditions in Ref. 8, 



The equipotential gate voltage on the boundary of a 1 
/.tm-thick Si02 dielectric Vg and the accelerating volt- 
age applied to the detector Vd are first used as partial 
Dirichlet boundary conditions in the solution of Poisson's 
equation 



(1) 



where the space charge density p is zero in the dielectric 
but equal to the donor impurity density from full deple- 
tion (calculated from room-temperature resistivity of 4 
kficm in the Si to match our device characteristics) . The 
dielectric constant e is also spatially dependent (=3.5 in 
the Si02 and =12 in the Si transport channel). The fi- 
nite differences method solves for the spatially discretized 
electrostatic potential (/) by a simple inversion of the ma- 
trix representation of the linear operator in Eqn. [U where 
Neumann boundary conditions (zero field normal to the 
boundary) are used on the boundaries not adjacent to 
the injector, detector, or gate equipotentials. An exam- 
ple electrostatic potential map calculated in this way for 
Vd =10V and Vq =2V is shown in Fig. [T](b), where we 



have additionally included a conduction band disconti- 
nuity of 3.1 eV between Si and SiO2.0 

Transport of a single electron spin component s in 
this potential is then governed by the classical two- 
dimensional drift-diffusion-relaxation equation (valid for 
the relatively high temperatures to be simulated and in 
an effectively bulk regime where the mean-free-path is 
much shorter than the device lengthscale) 



— = VDVs - V[i;s] - s/t + G, 
at 



(2) 



where the spin generation source G is nonzero only at 
the injector. We can ignore the other two vector com- 
ponents of spin because the empirical spin precession re- 
sults in a perpendicular magnetic field can be reduced 
to a transit-time distribution for just one component via 
Fourier transform.!^ |^ We therefore only need to de- 
termine the drift velocity vector field v and the diffu- 
sion coefficient scalar field D. Because the bulk spin 
lifetime in Si is several hundred ns at the temperature 
range of interest's'] and the effects we wish to model in- 
dicate much stronger spin relaxation than what is already 
known about the bulk, we set the spin lifetime r to be 
finite only at either the top and bottom interfaces as ex- 
plained below. 

From the potential <j>, we calculate the electric field 
£ = — V0. This is then used to extrapolate a 
spatially-dependent drift velocity vector v from empirical 
electric-field- and temperature-dependent time-of-fiight 
measurements '20'|. The (charge) diffusion coefficient can 
then be calculated using Einstein's relation D — fikgT/q 
where fi = \v\/\£\ is charge mobility, fcs is the Boltz- 
mann constant, q is the fundamental electron charge and 
temperature T—60K is chosen to match Ref. 

Whereas we can be confident in our application of bulk 
data to our model of the interior points of our Si chan- 
nel, transport at the interfaces is however greatly mod- 
ified by enhanced momentum scattering which reduces 
the mobility fi. 21| To determine the extent to which this 
occurs in our devices, we performed Hall measurements 
in magnetic field B perpendicular to the device plane 
with contacts to the lateral transport region as shown in 
Fig. [IJa). In contrast to ohmic transport, here with an 
undoped channel and rectifying Schottky contacts at the 
injector and detector, the carrier density is not constant 
but rather determined by the injected particle current 
density J and the transport time i a,s n — Jj^, where L 
is the transport length. This gives a Hall voltage of 



VHaL 



^ B wLB 
Jw — = — = — 

n t 



(3) 



where w/L ^ 1 gives the ratio of device width to trans- 
port length and we use the linear response approxima- 
tion L/t w fiVo/L. Hence, with ballistic hot electron 
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FIG. 2: (a) Plan view of 10-/im-thick lateral undoped Si spin 
transport device with Hall probes; (b) Example charge mobil- 
ity directly measured with transport of ballistically injected 
electrons in a perpendicular magnetic field with Vd = WV 
at lOOK using Eqn. [3] and the data shown in the inset. The 
two sets of data shown here demonstrate the reproducibility 
and lack of any spurious hysteresis. This strong dependence 
on gate voltage motivates the interfacial mobility reduction 
incorporated into our finite differences model. 



injection into otherwise insulating lateral semiconductor 
channels as we have here due to low doping, low tem- 
perature, and rectifying Schottky barrier contacts. 
the Hall voltage is surprisingly independent of injection 
current to first order. Furthermore, in contrast to Hall 
measurements on ohmic systems which yield the carrier 
density, here the magnetic field dependence (slope vs B) 
gives the mobility /x directly. 

The example experimental mobility data from appli- 
cation of Eqn. [3] shown in Fig. El^b) give evidence for 
a strong momentum scattering enhancement upon appli- 
cation of gate voltage, with a reduction by a factor of 2 
between Vc =0V and 4V. To account for this strong ef- 
fect in our model, we need to lower the bulk drift velocity 
value in elements adjacent to the interfaces in the trans- 
port channel by a significant factor. However, because of 
the low doping density and relative lack of bandbending, 
our electrons are not confined entirely to the interface 
layer as in a two dimensional electron gas. Therefore, 
the mobility reduction measurement is an average over 
some electrons which never touch the interface and it is 
difficult to assign a precise value to the interfacial mo- 
bility reduction in the finite differences model based on 
experimental results. Measurements on similar (Si/Si02 
field-effect) devices indicate a large characteristic mobil- 
ity reduction of »ilO in inversion layers due to e.g. inter- 
face roughness (2 ll. [23|: since this is consistent with our 
experimental observations, we therefore use that value in 



subsequent model calculations. 

The drift-diffusion-relaxation equation in Eqn. [2] de- 
scribes the evolution of spin density s in time and space. 
To model the spin transit time in a straightforward way, 
we can simulate the sourceless evolution of an initial delta 
function spin distribution at the injector which can be 
done by a suitable technique explicit in time such as the 



Crank-Nicolson method. [2J, |25| This corresponds to an 



iteration in time with stepsize At (starting with a delta- 
function nonzero only at the injector) of 



= (I - — H) 



(I- 



-H)s" 



(4) 



where s" is the array of spin density values at time 
t = nAt. H is the block-tridiagonal matrix representa- 
tion of the drift, diffusion, and relaxation linear operators 
(with inhomogeneous v and D) in the finite differences 
scheme and I is the identity matrix. Elements at all in- 
terfaces are treated appropriate to the boundary condi- 
tions chosen; insulating boundaries are used (including 
at the Si/Si02 interface due to the conduction band off- 
set barrier) except at the detector where s = so that 
spin flows out of the computational region to constitute 
a signal current. 

Although explicit iteration in time as in Eqn. U] is 
a possible route to calculating the transit time, due 
to causality constraints At must typically be chosen 
very small for numerical stability. Directly calculat- 
ing relatively long timescale transport phenomena in 
this way with large matrices necessary to discretize 
two-dimensional space can therefore be computationally 
costly due to the large number of iterations required. 

Therefore, instead of using the Crank-Nicolson scheme 
in Eqn. H] (or other time-explicit alternatives), we trans- 
form Eqn. [5] into the frequency domain to obtain an 
equation for s(w) — J s{t)exp{—iujt)dt. This allows 
us to directly simulate the steady-state response with 
spin precession in a perpendicular field (at frequency 
Lu — gfisB/h, where g is the electron spin g-factor, fig is 
the Bohr magneton, and h is the reduced Planck con- 
stant). Within the finite differences scheme, we then 
must simply perform a matrix inversion to solve 



(iujl - H)s G 



(5) 



where G is a constant source at the injector for each 
chosen value of to. Then, the time-domain response to 
an injected pulse can be reconstructed via the Fourier 
transform of the integrated transport current signal at 
the detector, Js{ijj) = -D^. 

Using discretization Aa;=200nm and Ay=50nm, we 
have implemented the finite difference scheme outlined 
above on our spatial domain using sparse matrix opera- 
tions in MATLAB. The results discussed below are insen- 
sitive to the choice of these discretization parameters so 
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FIG. 3: (a) Example steady-state spatial spin distributions 
with constant spin source at injector "i" obtained with ma- 
trix inversion of Eqn. [5]at uu = for several gate voltages (Vg) 
with detector voltage Vd =10V; (b) Spin transit-time distri- 
butions upon arrival at detector "d" as a function of gate volt- 
age at Vd =10V with Si/Si02 interfacial spin lifetime t —1 
ns. Inset: Example spin precession simulation from which the 
spin current transport distributions are obtained via Fourier 
transform; here Vg'=2V. 



long as the meshing is detailed enough to prevent compli- 
cations from spurious and unphysical oscillations which 
are well-known in finite-difference methods. !26ij In Fig. 
|3^a), we show an example steady-state spatial distribu- 
tion for cli = at detector voltage bias Vd=10V and gate 
voltage biases Vc=0.5, 2.5, and 4.5V (both relative to the 
injector source potential). The injected spin is clearly at- 
tracted from the injector source at the upper left of each 
image to the electrostatically-gated Si/Si02 interface at 
the bottom of each image, and the electron distribution 
fills a larger lateral region along that interface as the gate 
voltage increases. Eventually, Vd draws the electrons 
laterally to the detector, where it first passes through a 
region close to the exposed air/Si interface at the top of 
each image before being absorbed at the detector contact 
edge. 

Fig. El^b) shows the spin current transit time distribu- 
tion with a small interfacial lifetime r = Ins (only at the 
Si/Si02 backgate interface) reconstructed from simulated 
spin precession measurements by varying w (as shown in 
the inset) for several gate voltages 0.5V^ < Vg < 4.5^. 
Clearly, as the injected electrons are attracted more to 
the backgate Si/Si02 interface with larger gate voltage. 
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FIG. 4: Interfacial spin lifetime dependence of transport dis- 
tribution properties as a function of gate voltage Va- In 
(a), total spin polarization, (b) mean spin transit time, and 
(c) transit time uncertainty (distribution standard deviation) . 
Black lines are simulations which include relaxation only at 
the Si/Si02 interface, whereas grey dashed lines include re- 
laxation only at the air/Si top surface. Lifetimes of «1 ns 
are necessary to reproduce the experimentally-observed de- 
creasing trends of all these parameters'^, in contrast to the 
known bulk lifetime in Si of several hundred ns at temperature 
of 60K.[3 



several changes occur: the total spin polarization (area 
under the curve normalized by the t = oo result to ac- 
count for spin conservation) decreases, the mean tran- 
sit time decreases, and the distribution width (transit 
time uncertainty leading to spin dephasing) decreases. 
All these trends mimic the empirical results from Ref. 

This gate- voltage-dependent behavior is quantified by 
the distribution trends at different chosen Si/Si02 spin 
lifetimes r =0.1-10 ns shown by the solid lines in Fig. 
IHa)-(c). It can clearly be seen that for interfacial spin 
relaxation r > Ins, the distribution parameters (spin po- 
larization, mean transit time, and standard deviation) 
follow trends which are inconsistent with previous exper- 
imental observations in Ref. ^: in particular, the latter 
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two properties shown in Fig. llJJb) and (c) increase with 
increasing gate voltage when the hfetime is too long. 

If spin relaxation is instead included only at the top 
air/Si surface, the correspondence to experiment is even 
worse. As shown by dashed lines in Fig. IHa)-(c), the 
polarization, mean transit time, and transit time uncer- 
tainty are all predicted to rise with applied gate voltage, 
regardless of the chosen interfacial spin lifetime, contra- 
dicting experiment. Therefore, our prior experimental 
results in Ref. ^ cannot be explained by top-surface spin 
relaxation at the air/Si interfoce, nor purely electrostatic 
effects. Adequate simulation requires the inclusion of 
dominant Si/Si02 interfacial spin relaxation at rates that 
are distinct in magnitude from the intrinsic mechanisms 
in the bulk which would otherwise lead to a spin lifetime 
orders of magnitude longer than the transit times seen 
here.fi"] 

It should be noted that while this calculation is in- 
tended to capture the dominant transport mechanisms 
in thick-channel gated lateral devices, it is not an ex- 
plicit simulation of the exact experiment in Ref. @. For 
instance, the injector-to-detector distance in the model is 
smaller to keep the calculation tractable and we do not 
attempt here to include the effect of spin injection only at 
the interior of the source contact as in true ballistic hot 
electron injection with a tunnel junction cathode. How- 
ever, the features of the geometry are quite general, so 
modifications to account for these discrepancies are not 
expected to change the overall trends shown here. In ad- 
dition, due to experimental unknowns such as interface 
dipole formation due to gate doping level and type, etc, 
we do not have the ability to precisely determine the ex- 
act empirical electrostatic boundary conditions from the 
experimental voltages applied, as evidenced by our abil- 
ity to operate the device at negative gate voltage beyond 
where channel pinch-off conditions are expected. Q Our 
calculations here do however establish a consistent trend 
with increasing gate voltage which can only be explained 
by inclusion of a strong interfacial spin relaxation that 
becomes more and more dominant as the injected elec- 
trons are attracted to the Si/Si02 interface. 

Although we cannot yet directly identify the mech- 
anism of spin relaxation at this Si/Si02 interface, the 
strong reduction in charge mobility as a function of gate 
bias shown in Fig. suggests a contribution from the 
Elliott mechanism. [27[ In this case, because the two de- 
generate conduction band states are not pure spin states 
due to spin-orbit interaction, 0, 0] momentum scatter- 
ing - even from spin-independent potentials - can lead 
to spin flips. Closer proximity to the gate interface in- 
duced by rising gate voltage and its reduction in charge 
mobility from momentum scattering is therefore consis- 
tent with the lower interfacial spin lifetime we deduce 
from spin transport modeling. However, because the El- 
liott spin-flip rate is linearly related to the momentum 
scattering rate and hence inversely proportional to the 



mobility, a factor of two reduction in mobility over the 
gate bias range used here indicates this mechanism can- 
not be the dominant effect responsible for our observa- 
tions of spin lifetime suppression of two orders of magni- 
tude. The role of other relaxation mechanisms such as ex- 
change coupling to extrinsic paramagnetic spins (e.g. Pfc 



centers 



12l |) and electric field fluctuations causing spin- 
orbit effects [3] must be explored. Although recently 
the effects of transient spin-polarized electron trapping 
in bulk defects have been observed. [isjl the reduction in 
transit time with increasing bias seen in Ref. d is in- 
consistent with this mechanism playing a role here. It 
is hoped that further experimental efforts with e.g. spin 
resonance techniques on samples having deliberately pre- 
pared interfaces with variable defect density will be capa- 
ble of clarifying the details of such relaxation pathways. 
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